Stochastic Time-Dependent Current-Density-Functional Theory 



O 

o 

(N 



Massimiliano Di VentrgQ and Roberto D'AgostsQ 
University of California - San Diego, La Jolla, California 92093 
(Dated: February 6, 2008) 

A time-dependent current-density-functional theory for many-particle systems in interaction with 
arbitrary external baths is developed. We prove that, given the initial quantum state \^o) and 
the particle-bath interaction operator, two external vector potentials A(r,t) and A'(r,t) that pro- 
duce the same ensemble-averaged current density, j(r,t), must necessarily coincide up to a gauge 
transformation. This result greatly expands the applicability of time-dependent density-functional 
theory to open quantum systems, and allows for first-principles calculations of many-particle time 
evolution beyond Hamiltonian dynamics. 
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Time-dependent density-functional theory 

(TDDFT) [H is becoming an important tool in the 
study of the dynamics of many-particle systems [2^ . 
The theory, which extends the applicability of ground- 
state density- functional theory [2, 3] to non-equilibrium 
problems, rests solidly on a theorem proved by Runge 
and Gross (RG) The RG theorem shows that, 

for a fixed initial quantum state the external 

time-dependent scalar potential Vexti^ii) acting on the 
many-particle system is uniquely determined (up to a 
constant) by the time-dependent single particle density 
ri(r,t).^ It was subsequently realized by Vignale and 
Kohn 4] that, since the exchange-correlation (xc) scalar 
potential in TDDFT is a strongly non-local functional of 
the density, a theory formulated in terms of the current 
density j(r, i) allows for a local gradient expansion of 
the xc vector potential A.xc{i^,t). This theory goes under 
the name of time-dependent current density-functional 
theory (TDCDFT) 0, and has been applied to a 
variety of contexts ranging from optical spectra of solids 
and atoms , to dielectric properties of polymers [1] , 
to transport in nanoscale junctions 0, lloL fill , [l^ . to 
the study of systems with memory and/or dissipation 

Despite the enormous success that the above time- 
dependent density-functional based methods have been 
experiencing, one needs to acknowledge that they can 
strictly deal with particle systems evolving under Hamil- 
tonian dynamics. There is, however, a large class of phys- 
ical problems where one needs to consider interactions 
with an external bath. Examples include dephasing of a 
quantum system coupled to a boson field, spontaneous 
emission due to coupling with the zero-point energy fluc- 
tuations, non-radiative decay, etc. A TDCDFT able to 
deal with such open-system problems, without reverting 
to consider the microscopic dynamics of the bath, would 
be a tremendous asset in the study of many-particle prob- 
lems beyond Hamiltonian dynamics [2^. 

In this Letter, we introduce such a theory which 
we name "Stochastic TDCDFT" and prove a funda- 
mental existence theorem comparable to the theorems 
of TDDFT and TDCDFT. Our theorem allows us to 
rigourously map the many-particle dynamics in the pres- 



ence of an external bath into an effective single-particle 
dynamics in the presence of the same bath. 

In non-equilibrium quantum statistical mechanics, the 
effect of a bath is effectively described using stochastic 
processes: the fluctuations introduced by the presence of 
the bath are accompanied by dissipative effects in the 
particle dynamics. Our starting point is therefore the 
stochastic time-dependent Schrodinger equation, which 
reads (?i = 1) 



dt-^{v,t) = -iH{t)-^{r,t) ~ ]^V^V-^{r,t) 



Vt{t)^{v, t) 
(1) 

where € describes a stochastic process, and V is an op- 
erator which describes the bath, and its interaction with 
the many-particle system [26| . In the following we will 
assume that V is independent of time and uniform in 
space. These restrictions can be easily removed (pro- 
vided one adds the requirement that V admits a series 
expansion in time) and our theorem remains valid, but we 
maintain them for the sake of simplicity. A few examples 
of the operator V can be found in [l^ • The Hamiltonian 
operator is 



[K + eA(f,,t)]^ 
2m 



fi-fl), (2) 



with A{fi,t) a vector potential and U{fi—fj) a po- 
tential describing two-particle interactions With- 
out loss of generality the stochastic process, £{t), is cho- 
sen such that it has both zero ensemble average and 
(5— autocorrelation, i.e.. 



£(i) = 0; £{t)i{t') = 5{t~t'), 



(3) 



where the symbol indicates the average over a sta- 
tistical ensemble of identical systems all prepared in the 
same initial quantum state \'^o) f28j . 

The last term on the rhs of Eq. ([1]) describes precisely 
the "fluctuation" induced by the bath; the second term 
is the compensating "dissipative" part. Equation ([T]), 
with the stochastic process defined in (l3|), is a stochas- 
tic differential equation which preserves the ensemble- 
averaged norm and the orthogonality between any two 
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states, il^a{r,t) and ippiv.t), J dr 7/)* (r, t)V'/3(r, = <5a/3- 
Ensemble-averaged orthonormality can be easily demon- 
strated by integrating Eq. ([1]) in a small interval of 
time At, taking the scalar product between two states, 
and ensemble-averaging this product using the proper- 
ties © 0. 

Given any operator in the Heisenberg representa- 
tion, its average over the initial state is obtained via 
the standard quantum mechanical definition {0{t)) = 
{'9o\0{t)\^o)- If O is the density operator n, one can 
easily show that the dynamics induced by Eq. ^ pre- 
serves the ensemble-averaged particle number, i.e., the 
quantity N = J dr {n(r,t)) = f dr {^'^ {r,t)^{r,t)) is a 
constant of motion [16|. Physically, this corresponds to 
the case in which the bath may exchange particles with 
the system, while keeping their average number constant. 
Similarly, one can prove that, given any observable 0{t), 

its ensemble average time evolution, 0{t), satisfies [l6| 



dtO = 



H,0 



--v^vd--dv^v + v^dv. (4) 



In passing, we note that if [O, V^] = at any time, then 

the equation of motion ^ reduces to dfO — i H,0 

i.e., the presence of the bath does not affect, on aver- 
age, the evolution of the observable O. Since, in general, 
[H, V] ^ 0, the total energy is not conserved. Addition- 
ally, if p is the density matrix operator, from ([1]) we get 
the well-known (Lindblad) quantum master equation 



dtp = -i 



- Iv^Vp-^pV^V + VpVl 



We define the ensemble-average density 



n{r,t) = {n{r,t)), 



(5) 



(6) 



reasonable physical assumptions, given an initial condi- 
tion IvE'o), and an assigned bath operator another ex- 
ternal potential A'{r,t) which gives the same ensemble- 
averaged current density, must necessarily coincide, up 
to a gauge transformation, with A{r,t). 

Proof: To prove this statement we follow a line of rea- 
soning commonly used to establish similar theorems in 
TDDFT and TDCDFT [13, [II- Let us assume that the 
same ensemble-averaged density n(r, t) and current den- 
sity j(r,t) are also obtained from another many-particle 
system with Hamiltonian 



H'it) = y: 



2m 



lY^U'ih-f,), (9) 



evolving from an initial state |^'g) and following the 
stochastic Schrodinger equation ([1]) with the same bath 
operator V. |\E'q) gives, in the primed system, the same 
initial current and particle densities as in the unprimed 
system. 

The core of the demonstration is as follows: By writing 
the equations of motion for j(r,t) determined by both 
A(r, t) and A'(r, t), we obtain an equation of motion for 
the potential difference AA(r, t) = A(r, t) — A'(r, t). We 
then prove that AA(r, t) is completely determined by 
the initial condition via a series expansion in time about 
t = 0. Finally, if the two systems coincide then the unique 
solution is A A = 0, up to a gauge transformation. 

The equation of motion for the ensemble-averaged cur- 
rent density is easily obtained from the equation of mo- 
tion of the current density operator. From Eq. ^ we 
get 



9J(r,t) = !^dtA{r,t) ~ x [V x A{r,t)] 



m 



m 



and current density 



(10) 



j(r,<)^0-(r,t)), (7) 

where the current operator is defined as j{r,t) — 
with 



Pi + eA{ri,t) 



m 



(8) 



the velocity operator of particle i, and the symbol 
{A, B} = {AB -\- BA) is the anticommutator of any two 
operators A and B. 

We have now set the stage to formulate and demon- 
strate the following theorem. 

Theorem: Consider a many-particle system described 
by the dynamics in Eq. ([1]) with the many-body Hamil- 
tonian given by Eq. ([2]). Let n(r,t) and j(r,<) be 
the ensemble-averaged single-particle density and cur- 
rent density, respectively, with dynamics determined by 
the external vector potential A(r,t) and bath V. Under 



where we have defined 

g{r,t) = V^j{r,t)V - ij(r,t)yty _ iytf 
T{r, i) = - ^ (5(r - f,)VjC/ (r, - fj) + mV • a (r, t) 

(11) 



with the stress tensor cr(r, t) given by 



The first two terms on the rhs of Eq. pOj) describe the ef- 
fect of the applied electromagnetic field on the dynamics 
of the many-particle system; the third is due to particle- 
particle interactions while the last one is the "force" den- 
sity exerted by the bath on the system. 
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Equations similar to Eqs. ([8]) - ([12)) can now be written 
for the system with the vector potential A'(r,t). Simi- 
lar force terms T' and Q' appear in these new equations. 
T' and Q' differ from the same forces in the unprimed 
system, since the initial state, the external vector poten- 
tials and the velocity v are different. By assumption, the 
ensemble-averaged current and particle densities are the 
same in the two systems, thus 



atj(r,i) = !^9,A'(r,i) - x [V x A'(r,t)] 



TO 



(^'(r,t)) 



(0'(r,t)). 



(13) 



Taking the difference of Eqs. (fTO| and (fT3|) we arrive at 



n(r, t)5t Ayl(r, t) = j(r, t) x [V x AA(r, t)\ + AQ(r, t) 

(14) 

where AA(r,i) = A'(r,i) - A(r,t) and A(5(r,i) = 

Q'(r,i) - g(r,t) with Q(r,t) = (^(r,i)) + m(^(r,t)), 
and (5'(r, t) the same quantity but in the primed system. 

We now need to prove that Eq. p4| admits only one 
solution, i.e., AA(r, i) is completely determined by the 
averaged dynamics of the current and particle densities, 
once the coupling with the bath, V ^ is assigned. To this 
end we expand Eq. in series about t = and obtain 
an equation for the Z-th derivative of the vector potential 
AA(r, i). That one can expand this equation in a time 
series about < = follows immediately from the analyt- 
icity of the vector potential and Eq. ([4]) [2^. We thus 
arrive at the equation 

i-i 

«o(r)(Z + l)A^,+i(r) = - ^(fc + l)n,_fe(r) AA^+i (r) 

A;=0 

+ AQz(r) 

+ ^j,_fc(r) X [V X AAfe(r)] 

(15) 



a7(r,t) 



t=0 

does 



where, given an arbitrary function of time /(r, i), we 

have defined the series expansion fi{r) i 

We are now left to prove that the rhs of Eq. ([1 
not contain any term AAi+i(r). This follows from the 
fact that the dynamics of any ensemble-averaged oper- 
ator is given by Eq. (j4]). Indeed, this implies that the 
l-th time derivative of any operator can be expressed in 
terms of its derivatives of order k < I, time derivatives 
of the Hamiltonian of order k < I, and powers of the 
operators V and . The time derivatives of the Hamil- 
tonian do contain time derivatives of the vector potential 
A(r,t), but always of order k < I. Then on the rhs of 
Eq. (fT5|) no time derivative of order I + 1 appears. Equa- 
tion psp can be thus viewed as a recursive relation for 
the time derivatives of the vector potential AA(r, t). To 



complete the recursion procedure we only need to assign 
the initial value AAo(r) = A(r, i = 0) - A'(r,t = 0). 
Since in the unprimed and primed systems the densi- 
ties and current densities are, by hypothesis, equal, the 



initial condition is simply given by n(r, t 
(^o|jp(r-,t 



O)AAo(r) 



0)1^^0) ^{%\jpir,t = 0)1*;,), where jp(r) = 
(1/2to) X^ifei S{r—'ri)} is the paramagnetic current den- 
sity operator. 

The same considerations as in Ref. [3l about the finite- 
ness of the convergence radius of the time series (ITSt ap- 
ply to our case as well. We rule out the case of a vanishing 
convergence radius by observing that it seems implausi- 
ble that the smooth (in the ensemble-averaged sense) dy- 
namics induced by Eq. ([T]) can introduce a dramatic ex- 
plosion of the initial derivatives of AA(r). If this holds, 
the expansion procedure (|15p can be iterated from the 
convergence radius time onward. We have then proved 
that Eq. completely determines the vector potential 
AA(r, t) and thus, since A(r, t) is assumed known, it de- 
termines A'(r, t) uniquely, up to a gauge transformation. 

To finalize our proof, we consider the case in which 
U = U' and l^-o) = |*o)- ^ this holds, AAo(r) = 0. 
Then the recursion relation admits the unique solution 
A A/ (r) = for any Z, and at any instant of time t we have 
A(r,t) = A'(r,t) (still up to a gauge transformation). 

Discussion: The theorem we proved provides solid 
grounds for the development of a Kohn-Sham (KS) 
scheme for the study of the dynamics of open quan- 
tum systems, i.e., many-particle systems coupled to an 
external bath. Indeed, we can choose as the primed 
system a non-interacting one, i.e., U' — 0. We thus 
find that any (ensemble-averaged) current density that 
is interacting A-representable is also non-interacting A- 
representable [l^[3|- (It is worth to point out that, in 
general, the current density is not w-representable, i.e., 
the mapping between current density and external scalar 
potential is not invertible [20|.) The dynamics of the 
many-particle system can then be mapped into the dy- 
namics of an "auxiliary" KS Slater determinant evolving 
according to 

^^^^'■^(r, t) = - tH^^^{t)^^'^(r, t) - ]^V^V-^'^^{v, t) 



where 



-|-F£(t)«'^'^(r,t), 



[Pi + eA{fi,t) + eA^c{fi,t)]'^ 
2m 



(16) 



+ VH{n,t), 



„ (17) 

with Aa;c[j(r, i), l^'o), V^] the xc vector potential, and 
VniVjt) the Hartree potential. In actual calculations, 
one can apply, as a starting point, available approxima- 
tions for Axe, 0, [i^] but more work in identifying 
the role of the bath on the correlations of the system is 
necessary. 
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We conclude by noting that the present theorem shows 
that given a bath and its interaction with the system, the 
ensemble-averaged expectation value of any observable 
can be written in terms of the ensemble-averaged current 
density. The most simple example is the average total 
current flowing in the system. Indeed, given an arbitrary 
surface S we can easily show that the ensemble-averaged 
total current Is{t) that flows across that surface is 

is(i)= / JM)-dS= / j(^«)(r,i)-dS = /f^)(t), 
Js Js 

(18) 

i.e., it is equal to the ensemble-averaged KS current 

Ig^^\t), where j(^'^)(r,i) is the (ensemble-averaged) 
sum of expectation values of the one-electron current- 
density operator in the populated KS states ipf^{rA). 
The result of Eq. extends the one proved in Ref. [23| 
for finite and isolated systems [s^] to the case in which 
the system is coupled to a bath. However, we note that 
the theorem we have proved in this paper is not limited to 
any specific boundary condition on the current density, 
and is thus valid for infinite systems as well. 

In summary, we have introduced a TDCDFT theory for 
many-particle systems in interaction with arbitrary ex- 
ternal baths. We have named it Stochastic TDCDFT. We 
have proved the uniqueness (up to a gauge transforma- 
tion) of the external vector potential, once the initial con- 
dition and the bath operator are assigned. This theory 
greatly expands the applicability of density-functional- 
based methods to systems evolving beyond Hamiltonian 
dynamics. We thus expect it will be of great value in the 
study of several physical problems, ranging from optics 
to transport. 

We thank K. Burke for useful discussions. This work 
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